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K> , ABSTRACT 

■ We present results from the 1996 MUSICOS (MUlti-SIte Continuous Spectroscopy) 

campaign on the T Tauri star SU Aurigae. We find a 2.7-d periodicity in the He i 
(587.6 nm) line and somewhat longer, less well-pronounced periodicities in the Balmer 
lines and in Na D. Our observations support the suggestion that the wind and infall 
signatures are out of phase on SU Aur. We present Doppler images of SU Aur that 
have been obtained from least-squares deconvolved profiles. Images taken about one 
rotation apart show only limited overlap, in particular at low latitudes. This is in part 
due to limitations in signal-to-noise, and in part due to line profile deformations that 
arise from short-lived and/or non-surface features. The agreement at high latitudes is 
better and suggests that at least some longer-lived features are present. The analysis of 
Stokes V profiles yields a marginal magnetic field detection during one of the phases. 

Key words: stars: pre-main-sequence - stars: imaging - stars: magnetic fields - stars: 
individual: SU Aur 



1 INTRODUCTION 



* Based on observations collected at the Canada-France Hawaii 
3.6m telescope, the McDonald 2.1m telescope, the La Palma 2.5m 
Isaak Newton telescope, the 1.93m telescope at the Observatoirc 
de Haute-Provence, the Xinglong 2.2m telescope and the Univer- 
sity of Vienna's automatic photometric telescopes. 



1.1 T Tauri stars 

One of the outstanding problems in the study of the later 
stages of star formation is the nature of the mechanism by 
which proto-stars lose angular momentum as they contract. 
T Tauri stars (TTS) are a good test bed for various accre- 



2 Y.C. Unruh, et al. 



tion theories, in particular classical TTS which show strong 
accretion signatures. 

TTS are usually subdivided into three groups. The first 
group, classical TTS (CTTS), are very active pre-main se- 
quence stars of about 1 Mq or less. Their spectra arc typ- 
ically characterised by excess UV and IR emission cis well 
as strong, low-ionization emission-line activity. This is su- 
perimposed onto a nearly "normal" late-type stellar absorp- 
tion spectrum that is thereby diluted. This dilution can be 
so strong that virtually no absorption lines arc visible in 
the most extreme CTTS. Most CTTS spectra can be ex- 
plained by a model in which viscous accretion (at a rate of 
10~* to lO^^Mej/r"^) from dusty circumstellar discs with 
radii of tens to hundreds of AU (Beichman et al. 1986; My- 
ers et al. 1987) drives the emission line activity (Kenyon & 
Hartmann 1987; Cohen, Emerson & Beichman 1989). The 
way the star and its disc interact, however, is still some- 
what controversial. CTTS usually show very strong and er- 
ratic light-curve variations with no clear periodicities. They 
are typically rather slowly rotating with rotation periods of 
about 5 to 10 days or so and may possess hot spots due to 
accretion shocks as well as cool starspots. 

The second group is called weak-lino TTS (WTTS). 
They are not accreting and show no evidence for a circum- 
stellar disk (in the visible wavebands) . WTTS often display 
higher rotation rates than their more active counterparts 
with Prot being of the order of 1 to 3 days (though sec 
also Stassun ot al. (1999) who found no significant differ- 
ence between WTTS and CTTS rotation periods for stars 
in the Orion Nebula). Most WTTS have fight curves that 
show clear periodicities and can be explained by invoking 
cool star spots. WTTS that show no disk signatures in the 
infrared are sometimes also called "naked" T Tauri stars. 

The third group consists of early- type TTS (ETTS), 
sometimes also called UXORs (after the prototype UX Ori) . 
These arc stars of spectral type KO or earlier (Eaton & 
Herbst 1995). Strictly speaking, the UXOR objects include 
TTS as well as Herbig Ae/Be stars. Their light curves are 
characterised by strong and erratic dimmings (quite often 
by about 1 magnitude and more in the V band). The most 
common explanation for these dimmings is the obscuring 
of the stellar disk by optically thick circumstellar gas (see 
e.g. Grinin & Tambovtseva 1995). The presence of optically 
thick circumstellar disks around ETTS, however, has not 
yet been well established. A number of ETTS show infall 
signatures, usually in the form of transient redshifted ab- 
sorption lines. This has prompted speculation that the in- 
falling gas is due to evaporating solid bodies (e.g. Grady 
et al. 1996; Grinin et al. 1994). ETTS do not appear to 
show any changes in the veiling continuum. 

To date, only very few T Tauri stars have been Doppler 
imaged. The most popular targets have been the two WTTS 
V410 Tau (Joncour 1992; Joncour, Bertout & Menard 1994; 
Strassmeier 1994; Hatzcs 1994; Rice & Strassmeier 1996) 
and HDE 283572 (Joncour, Bertout & Bouvicr 1994; Strass- 
meier & Rice 1998). The two more active TTS that have 
been imaged are DF Tau (Unruh, Collier Cameron & Guen- 
ther 1998) and SU Aur (Petrov et al. 1996). 

Much effort has been invested recently in measuring the 
magnetic fields on TTS, mostly using Zeeman line broaden- 
ing. Typical magnetic-field measurements yield of the order 
of 1 to 2 kG for WTTS (see e.g. Basri, Marcy & Valenti 



Table 1. Parameters of SU Aur. 



Parameter 




Ref. 


spectral type 


G 2 


Cohen & Kuhi (1979) 


brightness (V) 


9 mag 


Cohen & Kuhi (1979) 


Mass 


2.25 Mq 


Cohen & Kuhi (1979) 




1.9 M© 


Dewarf, Guinan Sz 






Shaughnessy (1998) 


vsini 


« 60 km s-i 


Johns-Krull (1996), 






this paper 


Period (photom.) 


1.55, 2.73 d 


Herbst et al. (1987) 


Period (pliotom.) 


<3.5 d 


Bouvier et al. (1993) 


Period (Ha wing) 


3 d 


Jolms & Basri (1995b) 


Period (jJrad) 


3.03 ± 0.03 d 


Petrov et al. (1996) 


Period 


2.7 ±0.3 d 


this paper 



1992; Guenther & Emerson 1996). For CTTS, Johns-KruU, 
Valenti & Koresko (1999) determined a mean surface mag- 
netic field strength of 2.6 ±0.3 kG for BP Tau and Guenther 
et al. (1999) found mean magnetic fields of 2.35 ± 0.15 kG 
and 1.1 ± 0.2 kG for the T Tau and Lk Ca 15 respectively. 
These field strengths are within the range of values needed 
for magnetic accretion models, (though see Safier (1998) 
for a critique of the magnctosphcric infall model). Using 
spectro-polarimetry, magnetic fields can be detected directly 
by their Zeeman signatures (in this case no accurate spec- 
tral line synthesis is necessary). Using this technique, Donati 
et al. (1997) detected complex magnetic fields on the sur- 
faces of WTTS V410 Tau and HD 283572. Johns-Krull et al. 
(1999) found the He i (587.6 nm) emission line of BP Tau to 
be circularly polarised and inferred magnetic field strengths 
of about 2.5 kG. 



1.2 SU Aur 

SU Aur is a relatively bright T Tauri star that has been 
variously classified as a classical T Tauri star (see e.g. Gi- 
ampapa et al. 1993), or early- type T Tauri star (Herbst et al. 
1994). It has also been used as a prototype to define an 
SU Aur class of T Tauri stars (Herbig & Bell 1988). Some 
of its parameters are listed in Tab. 1. According to Ake- 
son et al. (2002), SU Aur has a rather low disk mass of 
log(Mi3)/M0) = — 5.1^g g. This is considerably lower than 
typical classical TTS disk masses. 

On account of its relative brightness, SU Aur has been 
the subject of several detailed spectroscopic studies (Gi- 
ampapa et al. 1993; Johns & Basri 1995b; Petrov et al. 1996). 
Johns & Basri (1995b) put forward a magnetospheric accre- 
tion model derived from the currently favoured model for 
CTTS. In the CTTS accretion model the disk is disrupted 
by the magnetic field of the star and accreting matter is 
channelled onto the star along magnetic field lines (Konigl 
1991; Cameron & Campbell 1993; Shu et al. 1994; Armitage 
& Clarke 1996). In the case of SU Aur, Johns & Basri 
(1995b) suggest that the magnetic dipole is slightly inclined 
with respect to the rotation axis of the star. Material from 
the disk accretes along magnetic field lines, but because of 
the inclination of the dipole, accretion is seen preferentially 
for one half of the stellar rotation, whereas a wind is ob- 
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Table 2. List of the observatories with the wavelength coverage 
and resolution of the spectrographs used. 



Observatory 


detector 


AA 


Resolving 






[nm] 


Power 


BAO {2.2m) 


1024^ Tektronix 


557 - 850 


45 000 


OHP (1.9m) 


1024^ Tektronix 


389 - 682 


45 000 


INT (2.5m) 


1024^ Tektronix 


486 - 849 


35 000 


MDO (2.1m) 


1200 X 400 Reticon 


547 - 673 


55 000 


CFHT (3.5m) 


2046^ STIS2 


410 - 814 


35 000 



served during the other half. This model has been dubbed 
the "eggbeater". 

Even though SU Aur has been monitored frequently, 
the photometric period determinations are very uncertain. 
This is because the light curve of SU Aur does not show 
any clear periodic variations (see also Sec. 2.3). To date, the 
most conclusive period determinations have not come from 
photometric measurements, but were derived from spectral 
variations in the emission linos (sec Johns & Basri 1995b, 
Petrov et al. 1996 and Tab. 1). In fact, most of our cur- 
rent knowledge stems from time-series investigations of the 
Balmer profiles, though, so far, observations always suffered 
from relatively large phase gaps. 

Despite the uncertainties in its period, SU Aur is, com- 
pared to most other classical and early- type TTS, a promis- 
ing candidate for surface-mapping by virtue of its small veil- 
ing and high rotation velocity. In this paper we present the 
observations of SU Aur taken during the 1996 MUSICOS 
campaign and give an overview of the emission-line and pho- 
tospheric line-profile variability. In Oliveira et al. (2000b) 
the emission and upper-atmospheric lines were analysed in 
more detail. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 MUSICOS 

The MUlti-SIte Continuous Spectroscopy (MUSICOS) net- 
work was set up to investigate the time-dependent behaviour 
of stars with rotation periods that necessitate coordinated 
observations from several longitudes (Catala et al. 1993). 

SU Aur was chosen as one of the targets for the MU- 
SICOS 1996 campaign that took place in November of that 
year. A total of 126 echelle spectra were taken over 10 nights 
at five different sites. At the Beijing Astronomical Observa- 
tory (BAO) in Xinglong, China, the setup consisted of a 
2.2-m telescope with an echelle spectrograph and a 1024 by 
1024 Tektronix detector. The 1.9-m telescope at the Ob- 
servatoire de Haute Provence (OHP, France) was used with 
the Elodie spectrograph (Baranne et al. 1996). On the 2.1- 
m Isaac-Newton- Telescope (INT) in La Palma, Spain, the 
ESA MUSICOS spectrograph (Baudrand & Bohm 1992) 
was linked to the Cassegrain focus with a fibre-optic ca- 
ble. At McDonald Observatory (MDO, Texas), we used the 
2.1-m Otto Struve Telescope with the Sandiford Cassegrain 
echelle spectrometer (McCarthy et al. 1993). The observa- 
tions at the Canada-France-Hawaii Telescope (CFHT) in- 
cluded observations in Stokes V. They were taken through a 
Cassegrain polarimetric unit, fibre-linked to the MUSICOS 



Table 3. Journal of observations of SU Aur from BAO in Xin- 
glong, China. The columns list the identifier for the observatory, 
the day (in 1996 November) on which the spectra were taken, the 
UT mid time, modified Julian date (HJD— 2450400), the expo- 
sure time and the number of spectra taken per group. The last 
column lists the approximate S/N ratio per pixel in the Ha order 
for a single exposure. In order to obtain an estimate of the S/N 
ratio per resolution element, these numbers have to be multiplied 
by about 1.3. 



Obs 


Date 


UT 


MHJD 


Exp time 


Group 
size 


S/N near 
Ha 


BAO 


19 


15:20 


7.14 


1800 


2 


80 


BAO 


19 


18:09 


7.26 


1800 


1 


80 


BAO 


19 


19:52 


7.33 


1800 


1 


85 


BAO 


20 


14:25 


8.10 


1800 


1 


50 


BAO 


21 


13:37 


9.07 


1800 


1 


50 


BAO 


21 


15:02 


9.13 


3600 


1 


75 


BAO 


21 


18:27 


9.27 


1800 


1 


55 


BAO 


21 


19:39 


9.32 


1800 


2 


60 


BAO 


24 


13:28 


12.06 


1800 


1 


80 


BAO 


24 


14:36 


12.11 


1800 


2 


80 


BAO 


24 


17:58 


12.25 


1800 


2 


75 


BAO 


24 


21:27 


12.39 


1800 


3 


50 



Table 4. Journal of observations for SU Aur from the OHP, 
Prance. In order to obtain an estimate of the S/N ratio per res- 
olution element in the Ha order, the numbers listed in the last 
column have to be multiplied by 1.5. 



Obs 


Date 


UT 


MHJD 


Exp time 


Group 


S/N near 












size 


Ha 


OHP 


18 


22:14 


6.43 


1800 


3 


25 


OHP 


19 


03:22 


6.64 


1800 


3 


35 


OHP 


20 


21:55 


8.41 


1800 


3 


40 


OHP 


21 


01:32 


8.56 


1800 


3 


55 


OHP 


21 


04:24 


8.68 


1800 


2 


50 


OHP 


23 


22:12 


11.43 


1800 


3 


45 


OHP 


24 


01:32 


11.56 


1800 


2 


50 


OHP 


24 


05:00 


11.71 


1800 


2 


30 


OHP 


24 


21:31 


12.40 


1800 


2 


50 



echelle spectrograph (Donati et al. 1999). The wavelength 
coverage and resolution of each telescope/spectrograph com- 
bination are shown in Tab. 2 (see also Catala et al. 1999). 

The journals of observations are given for each obser- 
vatory individually in Tabs 3 to 7. The exposure time listed 
is the time of each individual exposure, rather than that 
of the complete group. Similarly, the signal-to-noise ratios 
listed are average ratios per pixel for a single exposure in 
a group. The signal-to-noise ratio per resolution element is 
obtained after multiplication with the factor given in the 
caption of each table. 

The weather during the run was patchy, particularly at 
the CFHT, but we obtained good phase coverage between 
HJD 2450406.5 and HJD 2450410.0, and reasonable phase 
coverage for the subsequent days up to HJD 2450414.0 (see 
Fig. 1). 
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Table 5. Journal of observations for SU Aur from the INT, La 
Palma. The multiplication factor to obtain the S/N for each res- 
olution element is 1.4. 



Obs 


Date 


UT 


MHJD 


Exp time 


Group 


S/N near 












size 


Ha 


INT 


17 


02:08 


4.59 


1200 


1 


60 


INT 


19 


23:05 


7.46 


1800 


2 


55 


INT 


20 


03:56 


7.66 


1800 


2 


60 


INT 


20 


05:55 


7.79 


1800 


2 


70 


INT 


20 


22:40 


8.44 


1800 


2 


60 


INT 


21 


01:09 


8.55 


1800 


2 


55 


INT 


21 


03:25 


8.64 


1800 


2 


75 


INT 


21 


05:39 


8.74 


1800 


1 


85 


INT 


21 


23:18 


9.47 


1800 


2 


65 


INT 


22 


01:27 


9.56 


1800 


2 


65 


INT 


22 


03:35 


9.65 


1800 


2 


60 


INT 


22 


06:10 


9.76 


1800 


2 


60 


INT 


23 


23:14 


11.47 


1800 


2 


40 


INT 


24 


01:08 


11.55 


1800 


1 


60 


INT 


24 


01:47 


11.57 


900 


2 


40 


INT 


24 


03:49 


11.66 


670 


1 


35 


INT 


24 


21:50 


12.41 


1200 


3 


45 


INT 


25 


00:20 


12.51 


1200 


3 


60 


INT 


25 


03:26 


12.64 


1200 


3 


65 


INT 


25 


22:14 


13.43 


900 


2 


35 



Table 7. Journal of observations for SU Aur from the Canada- 
Prance-Hawaii Telescope. Columns as in Tab. 3, the multiplication 



factor for the S/N 


is 1.5. 










Obs 


Date 


UT 


MHJD 


Exp time 


Group 


S/N near 












size 


Ha 


CFHT 


19 


09:39 


6.90 


600 


4 


55 


CFHT 


19 


12:40 


7.03 


600 


4 


50 


CFHT 


20 


08:14 


7.84 


600 


4 


60 


CFHT 


21 


14:18 


9.10 


600 


4 


30 


CFHT 


26 


08:55 


13.87 


600 


4 


30 


CFHT 


26 


11:30 


13.98 


600 


4 


50 


CFHT 


27 


10:35 


14.94 


600 


4 


50 


CFHT 


27 


11:34 


14.98 


600 


4 


45 


CFHT 


27 


14:37 


15.11 


600 


4 


35 



CFH 
MDO 
INT I 
OHP 
BAO 



II I I 

III III III 
I II Mil 

III 

III I III 



MM 



III 



406 



408 410 412 

HJD - 2450000. 



Figure 1. Phase coverage during the MUSICOS campaign 



Table 6. Journal of observations for SU Aur from the McDon- 
ald Observatory, Texas. Columns as in Tab. 3, the multiplication 
factor for the S/N is 1.4. 



Obs 


Date 


UT 


MHJD 


Exp time 


Group 


S/N ne 










N 


size 


Ha 


MDO 


19 


05:36 


6.73 


1800 


1 


110 


MDO 


19 


06:03 


6.75 


1200 


1 


95 


MDO 


19 


08:05 


6.84 


1300 


2 


105 


MDO 


19 


10:33 


6.94 


1500 


2 


115 


MDO 


19 


12:14 


7.01 


1500 


1 


100 


MDO 


20 


06:05 


7.75 


1500 


2 


95 


MDO 


20 


08:38 


7.86 


1600 


1 


80 


MDO 


20 


09:10 


7.88 


1800 


1 


80 


MDO 


20 


12:02 


8.00 


1800 


2 


85 


MDO 


21 


06:20 


8.76 


1800 


2 


95 


MDO 


21 


09:04 


8.88 


1600 


2 


115 


MDO 


21 


11:15 


8.97 


1600 


1 


90 


MDO 


21 


11:46 


8.99 


1800 


1 


80 


MDO 


22 


05:36 


9.73 


1600 


1 


80 


MDO 


22 


07:42 


9.82 


1500 


2 


100 


MDO 


22 


10:17 


9.93 


1800 


2 


100 


MDO 


23 


05:33 


10.73 


1800 


2 


85 


MDO 


23 


08:34 


10.86 


1800 


1 


110 


MDO 


23 


09:03 


10.88 


1600 


1 


110 


MDO 


25 


05:24 


12.73 


1800 


2 


80 


MDO 


25 


08:08 


12.84 


1800 


2 


80 


MDO 


25 


10:39 


12.94 


1800 


2 


80 


MDO 


26 


04:14 


13.68 


1800 


1 


85 



2.2 Spectroscopic data reduction 

The (lata were reduced following standard procedures. Af- 
ter bias subtraction, each exposure was flat-fielded. For the 
data taken at the INT, BAO and MDO the orders were ex- 
tracted using the optimal extraction algorithm (Horne 1986) 



as implemented in echomop (Mills 1994), an echelle data 

reduction package distributed by STARLINK. Due to inter- 
nal reflections in the MDO spcctrograpii wc iiad to reject 
two orders in the 610 nm wavelength region. In the case of 
the OHP, the data are automatically reduced at the obser- 
vatory with the INTER- TACOS procedure (Baranne et al. 
1996) and were used as such. 

The reduction of the polarimetric data from the CFHT 
was performed using a dedicated procedure for extracting 
Stokes V and I parameters developped as part of the "Es- 
prit" data reduction package by Donati et al. (1997). Each 
exposure is a combination of four subexposurcs that have to 
be taken in order to get rid of background and instrumen- 
tal polarisation. The data reduction is described in detail in 
Donati et al. (1997). 

We have not attempted to flux-calibrate our spectra 
as this is notoriously difficult for cchcllc observations with- 
out simultaneous low-resolution spectroscopy or photome- 
try. When we refer to "flux" we are therefore not considering 
absolute flux values but the flux above (or below) the nor- 
malised continuum. Finally, as the resolution is different for 
each telescope, wc have rc-binncd all the single-line spectra 
shown here to a uniform resolution of A/AA = 30 000, cor- 
responding to a velocity resolution of about 10 km s~^. The 
convolved profiles used for the Dopplcr imaging in Sec. 4.3 
are binned to a resolution of A/AA = 40 000. 

2.3 Photometry 

Unfortunately, there is no contemporaneous light curve with 
good phase coverage. SU Aur was monitored during the 
1996/1997 observing season with the Automatic Photomet- 
ric Telescope (APT) of the University of Vienna (Strass- 
meier et al. 1997; Strassmeier, Serkowitsch & Granzer 1999), 
but only a few data points could be secured during the spec- 
troscopic campaign. The light curve is shown in Fig. 2. Typ- 
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Figure 2. The photometric lightcurve of SU Aur taken with the University of Vienna's APT in Fairbanks, Arizona. The x-axes show 
modified HJD, i.e. H,JD-2450000. The top graph shows the V-band Hghteurve over the course of almost 90 days. The bottom graph on the 
left shows the enlarged lightcurve during the MUSICOS campaign, but also including the recovery phase of the obscuration event. The 
middle plot shows V — R (lower curve, open squares) and R — I (top curve, solid diamonds) during the campaign. The lower right-hand 
plot shows V — R vs. V during the 1996/1997 observing season. The filled symbols are for data taken between HJD 2450405 and HJD 
2450425. 



ically for SU Aur, we can detect no clear rotational modu- 
lation in the lightcurve. Instead, it shows strong and erratic 
dimmings during which the star gets fainter by typically 
about 0.5 mag in the V-band within 3 to 5 days. The recov- 
ery is usually slower and can take 10 days or so (see Dewarf, 
Guinan & Shaughnessy 1998 and Nadalin, Dewarf & Guinan 
2000 for a description of similar events). This behaviour is 
typical for SU Aur and prompted Herbst et al. (1994) to 
classify SU Aur as an UXOr variable. SU Aur's dimmings 
are not quite as dramatic as those of other UXOrs that easily 
exceed differences of one magnitude in the V band. 

One such dimming event was observed starting from 
about HJD 2450408.5 during our spectroscopic monitoring 
campaign when the star's V-band magnitude decreased by 
more than 0.5 over the course of about 3 days. Unfortu- 
nately, the weather at the APT was such that only one 
night between HJD 2450410 and 2450416 was clear. At HJD 
2450417, the star was again dimmer by 0.6 mag (in V) than 
"usual" , but due to our lack of data we can not say whether 
it recovered between HJD 2450410 and 2450417. In agree- 
ment with what has been observed for other UXOrs, SU Aur 
becomes redder as it gets fainter (see lower right-hand plot in 
Fig. 2). However, during the 1996/1997 observing season, it 
does not show the "colour reversal" typical for ETTS, where 
the stars become bluer as they enter strong brightness min- 
ima (see e.g. Eaton & Herbst 1995). This could be due to 
the fact that the dimmings that we observed were not very 
extreme. As it is the disk that is believed to be responsible 
for the bluer colours of the systems once the star itself has 



been sufficiently occulted, the lack of colour reversal could 
also be due to SU Aur's disk being unusually faint. 



3 PERIOD DETERMINATIONS 

Photometric period determinations are difficult as SU Aur's 
brightness varies erratically (see Fig. 2). Photometric period 
measurements range from as short as 1.7 d (Nadalin, Dewarf 
& Guinan 2000) and 1.55 d to 2.73 d (Herbst et al. 1987) up 
to 3.4 d (Bouvier et al. 1993). For SU Aur, spectral Une- 
proffie variations seem to offer a better chance of tracking 
the star's rotation, though one has to bear in mind that 
a line such as Hq is formed at various heights in the stel- 
lar atmosphere and under a plethora of different conditions. 
The identification of Balmer-line variations with the surface 
rotation of the star therefore remains somewhat uncertain. 

Several period measurements have been obtained from 
analyses of the Balmer lines and of the He D3 line, all of them 
indicating periods around 3 d (Giampapa et al. 1993; Johns 
& Basri 1995b; Petrov et al. 1996). Johns & Basri (1995b) 
furthermore found that the behaviour of the red wing of 
H/3 and the blue wings of Ha and H/? were anticorrelated. 
This was interpreted as the signature of a magnetic infall 
and wind, in what they called the "eggbeater" model (see 
Sec. 1.1 and Johns & Basri 1995b) which has so far proven 
remarkably successful at explaining a range of accretion and 
wind phenomena. 

The measurements referred to above were taken from 
single-site observations. This introduces problems as the 
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life times of the rotationally modulated features seem to 
be short, perhaps due to unsteady accretion. Indeed, Smith 
et al. (1999) concluded that, for SU Aurigae, the dominant 
variability time scale for Ha is of the order of or shorter than 
1 hour. 

In the following we discuss the line profile changes and 
associated periods for Ha, H/3, Na D and He I (587.6 nm). 
Starting from 2-D periodograms we sum the fluxes in those 
velocity bins that show strong power peaks and use three 
different methods to search for periods. Note that adding the 
fluxes degrades the frequency resolution somewhat, though 
we took care to only sum velocity bands with very similiar 
appearance on the original periodograms. The main results 
of our period search are presented in Tab. 8 and the period- 
search methods are outlined in the following section. 

3.1 Period-search methods 

The three methods used to search for periods are the Lomb- 
Scargle algorithm, the "clean" algorithm and a technique 
based on the minimisation of information entropy. The 
Lomb-Scargle algorithm is well suited to irregularely spaced 
data and is described in Horne & Baliunas (1986) and Press 
et al. (1992). When applied to the summed fluxes, the false- 
alarm probability (FAP) that gives an indication of the like- 
lihood that a given power peak is due to random noise was 
calculated using a Monte-Carlo bootstrap method. For this 
we ram 1000 trials where we randomly selected fluxes from 
our original data set while leaving the number of data points 
and location (i.e. observing time) unchanged. We then cal- 
culated the periodogram for each synthetic data set, noting 
the maximum power for each trial. We plotted the probabil- 
ity for a given power value (in the power ranges obtained in 
the simulations) and extrapolated this to predict the proba- 
bility of the actual power at the period found for the original 
data set. In each case the highest values of the power found 
for the 1000 synthetic data sets were less than half of the 
power peaks in the observed data. 

Once a period has been determined, one obviously seeks 
to obtain an error estimate for this period. Most formulae for 
calculating the errors in the frequency (such as e.g. Mont- 
gomery & O'Donoghue 1999; Horne & Baliunas 1986) as- 
sume that the signal to be subtracted is sinusoidal. It turns 
out that this is not a good assumption for this data set. 
We therefore adopt a very conservative estimate and list the 
range of the periods that lie within the FWHM of the power 
peaks. The nominal error on the period is more than one 
order of magnitude lower than this. 

The second method that we used is the 'clean" algo- 
rithm (Roberts, Lehr & Dreher 1987) as implemented in 
the programme PERIOD (Dhillon & Privett 1997) distributed 
through STARLINK. The periodograms axe shown in Fig. 13 
and the uncertainty range in the period is listed in column 6 
of Tab. 8. The implementation of the "clean algorithm" that 
we used also calculates the FAP using the bootstrap method 
outlined above, however, it only gives numerical values in ex- 
cess of about 10~®, setting all other values automatically to 
zero. In our case, all FAPs were nominally "zero" , i.e. below 
10~^, and are therefore not listed in Tab. 8. 

Prompted by the non-sinusoidal shape of the He i 
(587.6nm) equivalent- width variations (see below), we de- 
cided to explore a third method for period searches which 
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Figure 3. Stacked Ha profiles. For clarity, we have grouped to- 
gether profiles that have been observed within 4 hours of each 
other. The continuum flux has been normalised and the y-axis 
offsets of the profiles corresponds to the modified HJD of the 
exposure times, i.e. (HJD - 2450400). 



relaxes the fundamental assumption of sinusoidal varia- 
tions, namely phase dispersion minimisation (e.g. Stelling- 
werf 1978). For our particular purposes we use the informar 
tion entropy minimisation technique of Cincotta, Mendez & 
Nuiicz (1995). In this scheme, several periods are tried and 
for each one the data are phased up on the given period. 
The data are normalized in amplitude and phase from to 
1. The idea is to measure how much of this normalized area 
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Figure 4. Top, right-hand side: The normaUscd Lomb-Scarglc 
periodogram of the Ha profiles of SU Aur. The left-liand axis 
shows the frequency (in d~^), tiic rigiit-liand axis the period 
(in days). The contours are normalised powers of 18, 26, 32 and 
40. The lowest contours correspond to a false-alarm probability 
(FAP) of about 10"®. To the left of the periodogram, we have 
plotted the window function for Ha. It has been scaled so that 
the power at zero frequency is unity. Underneath the periodogram 
the normalised mean (solid line) and variance profile (dashed line) 
arc shown. The normalised variance profile has been multiplied 
by a factor of 4 and offset by 1. It was calculated according to 
- - (Johns & Basri 1995a), where / is the 

mean profile. 

is filled with data points. Random signals will fill the area 

with data points while periodic signals will define a curve. 
The degree of space filling is quantified by the information 
entropy and periodic signals produce a minimum in the en- 
tropy at the appropriate period. In this method, the user 
chooses the period range to search, the number of trial pe- 
riods to explore, and the number of squares to divide the 
total normalized area (phase x amplitude) into. 

Typically we searched periods ranging from 0.5 to 10 
days at 2000 evenly spaced points and we found that 9 
squares produced the best results. The period uncertainty 
was again taken from the width of the minima and the sig- 
nificance of the period estimated using a Monte Carlo sim- 
ulation of the data. The results from this method are listed 
in the last three columns of Tab. 8. 

3.2 Ha 

SU Aur shows only Ha in emission. The higher Balmer lines 
arc in absorption, though they show varying degrees of filling 
in. Hq is a relatively "typical" wind profile with a marked 
blue-shifted absorption component that tends to reach be- 
low the continuum. It is highly variable on a number of 
time scales as illustrated in Fig. 3. Especially noteable are 
the redward-moving features that appear at a velocity of 
about -100 km s"^ around HJD 2450408 and HJD 2450412 
(see Oliveira et al. (2000a) for a discussion) and the strong 
enhancement of the blue emission wing at approximately 
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Figure 5. The flux in Ha above the (normalised) continuum in 
three different velocity bins. The diamonds show the flux varia- 
tions in the blue emission peak where the 3-day period was found 
(—300 to —210 km s~^). The triangles show the variations of a bin 
spanning -180 to -130 km where the periodogram suggests 
a 5-day period. The squares show the flux variations in the red 
wing where the periodogram also suggests a 5-day period (170 to 
250 km s~^). For clarity, the triangles and diamonds have been 
offset by -0.7 and -1.0 respectively. 

HJD 2450409 when the star was at its faintest during the 
campaign (see Sec. 6). 

The main plot in Fig. 4 shows the Lomb-Scargle pe- 
riodogram (see Horne & Baliunas 1986; Press et al. 1992) 
for Ha. The false-alarm probability (FAP) has been calcu- 
lated for each velocity bin according to Press et al. (1992). 
The FAP gives an estimate of how likely it is that a peak 
of a given height is due to random noise, a low value in- 
dicates a significant signal. The FAP is below 10~^ for all 
contours drawn in Fig. 4. We note that the formula given by 
Press et al. (1992) can underestimate the FAP considerably, 
in particular when the data points show strong clumping. A 
more accurate method to estimate the FAP is to run a boot- 
strap Monte Carlo simulation, as we have done for selected 
velocity bins (see Tab. 8). 

Along the frequency a^xis of each periodogram we have 
plotted the window function (as obtained from a discrete 
Fourier transform, since it is not possible to use the Lomb- 
Scargle analysis for windowed data). Underneath the peri- 
odogram we have sketched the mean profile (solid line) and 
the normalised variance profile (dashed line). Two periodic- 
ities at around 3 d and 5 d are apparent for Ha. The one at 
approximately 3 d is the one previously identified by Johns 
& Basri (1995b) and Petrov et al. (1996). While Johns & 
Basri (1995b) found the 3-day period in a velocity bin rang- 
ing from about -250 to -100 km s~^, we find a period of 
approximately 3 days further towards the blue, in a range 
from -300 to -170 km s~^. This 3-day period seems due to 
a weakening and strengthening of the outer blue wing. At 
the same time, we find a 5-day periodicity centred around 
a velocity of -150 km s~^ (where Johns & Basri (1995b) 
found the strongest evidence for a 3-day period) and over 
for almost all of the red emission wing. 

We suggest that the 5-day period is spurious, as the 
shape of the window function allows for substantial leakage 
into the 5-day bin (see Fig. 4). Furthermore, we do not find 
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Figure 6. Stacked H,3 profiles. For clarity, we have grouped to- 
gether profiles in 4-hour bins. The continuum flux has been nor- 
malised and the y-axis oflfeets of the profiles correspond to the 
observing times (HJD - 2450400). The broadened profile of the 
Sun has been subtracted from all spectra. 



much evidence for a 5-day period when the intensities of the 
Ha profile are plotted against time in the regions where pe- 
riodicities are suggested in the periodogram (see Fig. 5). We 
found that exclusion of the profiles between HJD 2450409 
and HJD 2450410, that show the strong flux enhancement, 
leaves the position of the power peaks largely intact, so that 
it is unlikely that the flux enhancement alone is responsible 
for the 5-day period. 
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Figure 7. The normalised Lomb-Scargle periodogram of the H/9 
profiles of SU Aur (top figure, RHS). The y-axis show the fre- 
quency (on the left) and the period (on the right). Also shown is 
the window function for H/3 (top left-hand plot). The FAP of the 
lowest contour is 2 x 10~^, the FAP of the second contour (at a 
normalised power of 15) is 10~®. The bottom figure shows the H/3 
mean profile (solid line) along with the variance profile (dashed 
line). In both cases the rotationally broadened solar H/3 profile 
has been subtracted. The variance profile has also been expanded 
by a factor of 4 and oflfeet by 0.5. 



3.3 H/3 

Because of its different source function, H/3 is formed over 
a smaller range of temperatures and densities than Ha. 
Theoretical studios (sec e.g. Hartmann, Howctt & Calvet 
1994; Calvet & Hartmann 1992) as well as observations (Ed- 
wards et al. 1994), have shown that H/3 reveals the signature 
of mass infall much more readily than Ha. This is indeed also 
the case for SU Aur. 

Fig. 6 shows the stacked and grouped profiles of H/3 
during our run where we have subtracted the rotationally 
broadened H/3 proflle of the Sun. We were able to observe 
H/3 simultaneously from only three of the sites (CFHT, INT 
and OHP) and therefore obtained rather poor phase cover- 
age. The periodogram for H/3 is shown in Fig. 7. In agree- 
ment with Johns & Basri (1995b) we find that the period 
of approximately 3 days now occurs in the red as well as in 
the blue part of the profile. The periodicity in the red wing 
is not only apparent in the periodogram, but is, in contrast 
to Ha, now also visible in the stacked spectra (Fig. 6). As 
already suggested by the lower and less well-defined peats in 
the periodogram, the periodicity in the blue wing is harder 
to trace. 

Johns & Basri (1995b) suggested that the red and blue 
wings of the H/3 line are anti-correlated. Because of a lack of 
phase coverage and the presence of the small moving emis- 
sion bumps in the blue profile wings, it is hard to see this 
directly in Fig. 6, though more detailed analysis essentially 
confirms their original conjecture. We refer to Oliveira et al. 
(2000b) for plots showing the anti-correlation between the 
red and blue wings. 
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Figure 8. Top: The normalised Lomb-Scarglc pcriodogram of 
the Na D profiles of SU Aur (RHS) and the window function of 
the data set (LHS). The lowest contour corresponds to a FAP of 
7 X 10"®. Bottom; The mean and variance profiles for the Na D 
lines. 



3.4 Na D doublet 

Due to the interstellar absorption lines, light pollution and 
lower intrinsic variability, the sodium doublet is slightly 
more diflicult to analyse. However, they still show a num- 
ber of similarities with the Balmer lines, including redward- 
moving features in their blue absorption wings (see Oliveira 
et al. (2000a) for a detailed analysis). The pcriodogram for 
the Na D lines is shown in Fig. 8. We recover once again 
a period of slightly less than three days in the red parts 
(from about 30 to 130 km s~^) of the profile. In the same 
velocity region, we can now also see a strong secondary pe- 
riod at about 1.4 days. The variability at the line centre is 
strong, but does not show any periodicities. As with Ha, 
there is also a (probably spurious) power peak close to 5 
days, though now only in a narrow velocity range centered 
at -70 km s"^ 



3.5 He I 587.6 nm 

The stacked profiles of the He I 587.6 nm line are shown in 
Fig. 9. This plot clearly shows the dramatic changes in line 
strength and the mostly asymmetric shape of the line profile. 
Oliveira et al. (2000b) argued that the He line can most eas- 
ily be understood in terms of a superposition of two absorp- 
tion components, one at rest velocity, and one red-shifted by 
about 80 km s~^. Note that while the minimum observed 
He equivalent width of about 50 mA is not too dissimilar 
to what is usually observed on G dwarfs (Saar et al. 1997), 
the maxima of up to 600 mA are quite exceptional. Solar 
observations show that the He D3 absorption increases in 
strength with increasing activity, but then starts to fill and, 
for the largest flares, turns into emission (Svestka 1972). 
Bray (1964) reports equivalent widths of up to 400 mA in 
solar flares, but more recent observations seem to have fo- 
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Figure 9. The stacked He I profiles, grouped together in 4-hour 
bins. As the line-profile changes in the helium line are not very 
large, they have been amplified by a factor of 4 for this plot. 
(Despite other appeajrances, the noise level for the He i profiles is 
in fact very similar to the noise level for Na D). 



cussed on large and off-limb or near-limb flares that then to 
show He D3 in emission. 

In their Hc-line-formation calculations Andretta & Gi- 
ampapa (1995) identify two distinct formation regions in 
dwarf atmospheres, namely in the upper chromosphere and 
in the much hotter plateau region. In the density regimes 
considered by Andretta & Giampapa (1995) the "turnover" 
where the equivalent width starts to decrease due to reversal 
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Figure 10. Top: The normalised Lomb-Scargle periodogram of 
the He I profiles of SU Aur and the window function. Peaks higher 
than the lowest contour have a FAP of less than 4x10^**. Bottom: 
The mean and variance profile of He I over the run. The variance 
profile (dashed line) has been amplified by a fa<;tor of 4 and ofTset 
by 0.5. 



is never reached and they estimate a maximum equivalent 
width of 100 to 150 mA. This is expected to increase for sub- 
giants. Note that the strength of the He D3 hno is mainly 
electron density rather than temperature sensitive. 

The 2-D Lomb-Scargle periodogram for He D3 is shown 
in Fig. 10. It shows a very clear periodicity at 2.7 d with a 
strong secondary peak at 1.4 d for velocities ranging from 
about —30 to 160 km s~^. Pctrov ct al. (1996) observed 
periodic or quasi-periodic radial-velocity variations in the 
Balmer lines and in He D3. We find that there is indeed a 
trend for the radial velocity to be further redshifted as the 
equivalent width of the He line increases. But because of 
the line asymmetry, the velocity shift of the lino centre is 
not very well determined. The periodically varying equiva- 
lent widths of the absorption, however, are readily appar- 
ent, as illustrated in Fig. 11. The variations are clearly non- 
sinusoidal and display very strong and sharp peaks for the 
largest equivalent widths. This is what first prompted us to 
use the phase-minimisation technique which indeed places 
somewhat narrower constraints on the period of the He I 
variations (see Tab. 8). 

The period-dependence of the entropy is plotted for the 
He I line in Fig. 12 and shows a very pronounced minimum 
at 2.7 d. The secondary period of about 1.4 d is also recov- 
ered, though it is less prominent than in the corresponding 
"clean" cd or Lomb-Scargle poriodograms. The dashed line 
shows the mean entropy for 1000 independent Monte-Carlo 
trials. The 1000 individual trials look very similar to the 
mean shown in the figure. For each trial, the entropy mini- 
mum was recorded. The distribution of the 1000 trial min- 
ima peaks at an entropy of 0.8745 with a fairly sharp cutoff 
above this value and a Gaussian tail extending to lower en- 
tropy values. None of the trials produced a minimum lower 
than 0.843, while the original He I data shows a minimum 
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Figure 11. The He I equivalent width as a function of time. 

The different symbols show measurements from the different ob- 
servatories, namely the BAO (plus signs), OHP (crosses), INT 
(diamonds), MDO (triangles) and CFHT (squares). 
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Figure 12. The solid line shows the entropy as a function of 
period for the He i equivalent width (see text). Also shown is the 

mean entropy for 1000 Monte-Carlo trials where the time spacing 
has been kept while the order of the equivalent width data points 
has been randomized. 



entropy of 0.675 at a period of 2.7 days. From a Gaussian 
fit to the peak and lower values of the trial distribution we 
find that the 2.7-day period is significant at the > 21cr level. 

We find essentially the same periodicity in both compo- 
nents of the He i line, though the shape of the periodogram 
is somewhat different in that the variability is stronger and 
the secondary peak is more pronounced in the main absorp- 
tion line than in the redshifted component (see Fig. 13d). 
Note that Oliveira et al. (2000b) found that the variability 
of the redshifted absorption correlated well with the red- 
shifted absorption of the NaD lines, but lagged behind that 
of the helium rest- velocity component. 



3.6 Comparison of the different lines and methods 

The results from our period searches are summarised in 
Tab. 8 and, for the "clean" algorithm, are also plotted in 
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Fig. 13. While there are small differences in the periods 
found for the different hnes with the periods recovered for 
the Ha line being the largest (2.8 and 2.9 d), the widths of 
the power peaks are such that a single period is compatible 
with all our observations. Wc favour the slightly shorter pe- 
riods obtained in lines other than Ha, mainly because of the 
large range of regimes where Ha is formed, some of which 
may be well away from the stellar surface. There may hence 
not be sufficient magnetic "locking" to force co- rotation and 
a single wcll-dcfincd period. 

The widths of the power peaks do not change noticeably 
when isolated measurements such as the first INT exposure 
or the last three CFHT exposures are excluded. Wc suspect 
that this is due to several features of our data and the star, 
the main ones being that we have only covered two to three 
rotation periods and that the variations are generally non- 
sinusoidal. Note that the entropy minimisation technique 
that does not assume sinusoidal variations, shows a much 
narrower peak than either the Lomb-Scarglo or the "clean" 
method. Furthermore, we find that the window function al- 
lows for leakage of power at a frequency of about 0.37 d~^, or 
a period of 2.6 d, which is uncomfortably close to our main 
period. Inspection of the stacked profiles also shows that 
while profiles that arc separated by approximately one pe- 
riod do show similarities in the velocity range of the power 
peaks, they are certainly subject to strong variations be- 
tween one period and the next. We are hence dealing with 
a star with very strong intrinsic variability. One also has to 
keep in mind that the strong dimming of SU Aur may have 
introduced an additional trend in the fluxes. One of the char- 
acteristics of ETTS is the erratic dimming with simultaneous 
increases in the Ha relative flux. Unfortunately, our photo- 
metric time coverage is too sparse to allow cross-correlations 
between the photometric fluxes and the spectroscopic data. 



4 PHOTOSPHERIC LINES AND DOPPLER 
IMAGES 

4.1 Least-squares deconvolution 

One of the original aims of the MUSICOS campaign was 
to obtain Doppler images of SU Aur and to locate possi- 
ble infall sites for matter accreting from the stellar disk. 
The presence of spots on the stellar surface should produce 
line-profile deformations in the photosphcric lines. However, 
in our case, the signal-to-noise ratio is not good enough to 
allow us to detect these deformations in individual photo- 
spheric linos. One possibility is to combine the information 
contained in all the photospheric lines in the observed spec- 
tral range so as to increase the S/N. We have done this us- 
ing the technique of least-squares deconvolution (LSD, see 
e.g. Donati et al. 1997; Barnes et al. 1998). As we have gath- 
ered data at five different observatories with different setups, 
it was not possible to use the same lines in the deconvolution 
for each observatory. 

One implicit assumption when using LSD profiles for 
Doppler images is that all lines used in the deconvolution 
should show a similar response to surface inhomogcncitics. 
To test this, we generated LSD profiles with different line 
subsets, such as e.g. singly ionised lines, "light" elements, 
lines in the red or blue parts of the spectrum only or strong 
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Figure 13. Top to bottom: The "cleaned" pcriodograms for Ha, 
H/3, Na D and He. The periodogram for Ha shows the power 
over several velocity ranges, namely —300 to —210 km s^-' (solid 
line), -180 to -130 km s^^ (dotted line) and 170 to 250 km s"! 
(dashed line). These bins correspond to the ones chosen for the 
flux plots in Fig. 5 and also represent the bins with largest power 
in the 2-D Lomb-Scargle periodogram. The H/3 periodogram is 
for the average flux between 60 and 185 km s~^; taJiing a smaller 
bin width results in a very similar periodogram with the same 
main period. The Na D periodogram is for flux between 50 and 
70 km s~^. For He I, the solid line is for the average flux between 
20 and 50 km s~^, while the dashed line is for an average taken 
from 100 to 120 km s'^. 



linos only. This was done for the CFHT data only, as it 
provided a large spectral range and the best S/N. 

Subsctting for different elements or ionisation stages is 
not straightforward as neighbouring blends can easily skew 
the profiles. The only element lighter than calcium that 
showed relatively unblended lines in our wavelength range 
was silicon and we were only able to use 4 lines. The dashed 
line in Fig. 14a shows the profile for the 4 silicon lines for the 
CFHT. The solid line is the LSD profile for all elements heav- 
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Table 8. Table listing the periods found in the normalised Lomb-Scargle periodograms (column 3), from the flux variations using the 
"clean" algorithm (column 6) and the entropy minimisation technique (column 8) . The 2nd column gives the range of velocity bins that 
have been co-added for the periodogram. For the periods calculated using the Lomb-Scargle algorithm, we have listed the false-alarm 
probability (FAP) according to a Monte-Carlo bootstrap method. The FAP for the "cleaned" periods is below 10~® for all periods (and 
has not been listed, sec text). The ranges listed in columns 4 and 7 give a very conservative estimate on the period, assuming that the 
width of the power-peak at half maximum can be used as an indication of the possible spread in the period. For the entropy minimisation 
the error on the period as well as its signiflcance are listed in columns 9 and 10. The last row shows the period for the equivalent width 
variations of the He i line. 
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ier than and including calcium. Both profiles are "mean" 
profiles wliere the six exposures with tfie fiigliest S/N have 
been added. The mean silicon profile is very asymmetric 
with more absorption in the red part of the line. This is 
not due to the small number of lines that were used in the 
deconvolution or due to blending. In fact, all the strong sil- 
icon lines were inspected individually and showed the same 
asymmetric shape. Catala et al. (1999) found very marked 
difTerences between different ionisation species and elements 
for AB Aurigac, where an extra emission component was 
seen in the higher ionised species. Unfortunately, the singly- 
ionised lines in SU Aur are too weak and too blended in 
order to carry out a similar analysis. 

As far as the line variability and the line-profile defor- 
mations arc concerned, it appears that there is no difference 
between the various subsets or indeed the LSD profiles ob- 
tained using all available lines, so that using all lines in the 
spectrum should be a valid approach. Unfortunately, the 
S/N level of the smaller subsets is too low to quantify the 
agreement (or indeed the disagreement) between them. Wo 
decided to err on the side of caution and exclude possible 
artefacts due to the lighter elements and used only neutral 
lines of the heavier atoms for the Doppler images presented 
belowi . 

While least-squares deconvolution is necessary here to 
achieve a high enough signal-to-noise ratio, it introduces 
problems not usually encountered in multi-site Doppler 
imaging. This is due to the very different spectral ranges 
and the different setups at the 5 sites. This is shown in the 
middle graph of Fig. 14 where the mean profiles of all 5 sites 
are plotted. It turns out the the mean profiles at the INT and 
CFHT agree reasonably well, but that the profiles at BAO 
and MDO look very different. Given that not all phases are 
covered equally at the different sites, small differences in the 
mean profile can be expected. Nevertheless, the differences 
that we observe are too large to be attributed to a different 
phase coverage. 



We took care to exclude a large-enough section around the re- 
jected lines to avoid artefacts through blended lines. 



Table 9. The number of lines used at the different observatories. 
The first column indicates the observatory. The second column 
gives the number of different lines used for the convolution. The 
third column lists the number of line "images" that were actually 
used for the calculations. This number tends to be higher than 
the one given in column two as some lines were recorded twice on 
the chip. The last column gives the weighted central wavelength 
of the final profile. Note that we used the same line list for the 
CFHT and the INT. 



Observatory lines images Ac 

[nm] 



BAO 


352 


400 


656.7 


OHP 


1828 


2160 


526.7 


INT 


950 


1290 


584.8 


MDO 


311 


433 


601.4 


CFHT 


1728 


2660 


563.3 



We suggest that the differences are mostly due to the 
different setups used at each observatory for the following 
reasons. The number of lines registered for each spectro- 
graph setup is different at each observatory. Therefore, even 
using the same starting list of lines can produce different 
mean profiles. This is because some lines are registered twice 
at some spectrographs or fall into gaps in the spectral cov- 
erage. The number of lines included in the LSD line list is 
given in column 2 of Tab. 9 along with the actual number 
of line images (column 3) and the "central" wavelength of 
the deconvolved profile (column 4). For the deconvolution, 
each line is weighted according to its line depth, giving least 
weight to the weakest lines. Furthermore, lines on the edge of 
the chip whore the continuum flux level (and hence the S /N) 
is low are also given reduced weight. This latter weighting 
differs for each observatory. Wo tested the influence of this 
weighting for data from the CFHT and the INT. As the over- 
lap between the spectral ranges at these telescopes is very 
high, wc could use the same line list to produce two differ- 
ent mean proflles at each observatory, one of which takes the 
S/N into account. The mean profiles at the INT and CFHT 
that neglect the error weighting are indeed very similar so 
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Figure 14. Top figure: Plots of diflferent mean LSD profiles from 
the data taken at the CFHT. The solid line is the profile that has 
been obtained by deconvolving all neutral lines from elements 
heavier than (and including) calcium. The dashed line shows the 
mean profile of the four strongest unblended silicon lines. Mid- 
dle figure: The rather diff'erent mean profiles at the 5 observato- 
ries. Bottom figure: Comparison of LSD profiles taken at approx- 
imately the same time but at different observatories. The profile 
deformations agree reasonably well between the different obser- 
vatories, though there are some important discrepancies for some 
observations. The uppermost profiles were observed around modi- 
fied HJD {MHJD=HJD-2450400) of 12.4. The thick and thin dot- 
dashed profiles show two profiles from the MDO taken about 3 
hours apart, corresponding to a phase difference of only 0.04. This 
phase difference is too small to show the motion of the surface 
features. The difference between the two profiles is probably due 
to noise. The middle profiles were observed between MHJD 8.7 
and MHJD 8.8 and the bottom profiles between MHJD 7.7 and 
MHJD 7.9. 
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Figure 15. A stacked plot of the residuals of the deconvolved 
profiles. White is at a flux level of 1.008 and black at a flux level 
of 0.992. 



that the same template lookup table (see also Sec. 4.3) can 
be used for the INT and CFHT profiles. 

More worrying than the small shape difforonccs in the 
mean profiles of the different observatories arc occasional 
differences in profile deformations observed at roughly the 
same phase that can not be attributed to noisy data. This 
is illustrated in the bottom plot of Fig. 14, where profiles 
that were taken almost simultaneously at different sites are 



compared 



4.2 Photospheric line-profile variations 

We find marked shape changes in the photospheric line- 
profiles as mirrored in the convolved profiles. This is shown 
in the grey-scale plot of the residual deconvolved line profiles 



^ Note that the agrecincrit between individual strong lines, such 
as, e.g.. Ha that were observed at different observatories but at 
similar times is very good. This suggests that the discrepancy 
arises mainly because of the low S/N. 
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in Fig. 15. The residuals liave been produced by subtract- 
ing the mean profile of each observatory from the individual 
exposures at this observatory. Fig. 15 illustrates the mar 
jor problem for Doppler imaging SU Aur. Despite the clear 
periodicities found in the He I D3 and the Na D lines, no ob- 
vious periodicity appears to be present in the deconvolved 
photospheric lines. In fact, a periodogram analysis of the 
(mean-subtracted) deconvolved profiles does not show any 
significant peaks on the stellar rotation time scale. We do not 
think that this is due to the differences in the mean profiles 
at the different observatories, as a periodogram analysis of 
the CFHT and INT profiles that were produced with iden- 
tical weighting and that agree generally well did also not 
show any periodicities. This begs the question of whether 
large contrasting surface features are present on SU Aur or 
whether the line-profile changes are primarily due to causes 
other than star spots. 

One such cause could bo the strong dimming (usually 
attributed to circumstellar dust obscuration) that we ob- 
served around HJD 2450410 (see Fig. 2). The time of the 
faintest observation coincides reasonably well with one of 
the peaks in the He equivalent width and it is preceded by 
a very strong fiux increase in the Ha line. The behaviour of 
the LSD profiles during the dimming is somewhat ambigu- 
ous, though the observations at the MDO that cover this 
time span best, suggest a decrease of the equivalent width 
from 5% above average at HJD 2450409 (Ha maximum) to 
5% below average at HJD 2450410 (He maximum). Note 
that the scatter of the LSD equivalent widths is of the order 
of 3% at the MDO. If the variations in the light level were 
veiling-like, e.g., if they arose from excess emission due to hot 
spots, low light levels should entail increased photospheric 
line widths which does not fit in well with our observations. 

The predictions for the He line depend on what line- 
formation scenario is assumed. If the He absorption is pro- 
duced in the hot spot, the He equivalent width should be at 
a maximum during large veiling, though there might also 
be a dilution effect due to the excess emission. If, how- 
ever, the spot produces He emission and line-infilling, we 
would expect the He equivalent width to anticorrelate with 
the broad-band emission. While we lack high-time resolu- 
tion photospheric observations, a comparison between SU 
Aur's light curve and the variations in the He 587.6 nm line 
does suggest that their variability time scales axe different 
and that their variability is not correlated. This favours a 
scenario where the periodic He equivalent width variations 
are produced in dense infall columns or hot spots while the 
more erratic photometric variability is mainly due to dust. It 
is likely that the photospheric equivalent widths are due to 
a combination of both effects thus masking any clear trends 
and periodicities. 

We checked for indications of a colour change as a con- 
sequence of the photometric dimming. For the INT which 
offers a relatively good time and wavelength coverage, there 
is no discernable trend distinguishing the lines redward and 
blueward of 650 nm. The high-S/N data of the CFHT did 
also not show any colour trend. This is in good agreement 
with the findings by Herbst & Shevchenko (1999). While we 
do not see any colour changes during the dimming event, 
we do find that as SU Aur became fainter, i.e. from around 
HJD 2450408.0 onwards, all profiles show a strong asym- 
metry with the flux in the red part of the proflle being 



P = 2.7 d P = 2.7 d 





velocity [km/s] velocity [km/s] 

(a) (b) 

Figure 16. The left-hand figure shows the residuals as in the 

previous figure, but binned according to a period of 2.7 d. The 
right-hand figure is similar, except tliat only the profiles from the 
CFHT and INT have been included. We have also excluded the 
first profile at the INT (taken at HJD 2450404.6, well before the 
remainder of the data) and the profiles during the "dimming" . 
The grey scale runs from a fiux level of 0.992 (black) to 1.008 
(white) for the left-hand figure and from 0.995 to 1.005 for the 
right-hand figure. It is evident that some profiles do not fit in well 
with adjacent profiles. This is in part due to noise (see also the 
previous figure) , but also due to a lack of repeated structure and 
the insecurity in the stellar period. 



suppressed. It appears that this effect is strongest for the 
singly ionised linos while hardly visible for neutral lines with 
low atomic numbers (such as sodium, magnesium and sili- 
con). The profiles remain very similar for at least 1.2 d (see 
Fig. 15). As it takes roughly one third of a rotation period 
for a spot to cross the visible hemisphere of the star, this line 
asymmetry is clearly too long-lived and too stationary to be 
produced by a surface spot. Whether the asymmetry and 
the dimming are really linked is difficult to establish as the 
photometric observations are only very sparcely sampled. 

Whatever the cause of the profile asymmetry, recon- 
structions based on the complete data set will run into diffi- 
culties. We therefore reconstructed images where the profiles 
between HJD 2450407.75 and 2450409.0 have been excluded 
(sec next section). Even then we struggle with an appar- 
ent lack of periodicity. This is visible in the phase-binned 
residuals shown in Fig. 16, where some of the neighbouring 
spectra show very different features. Any selected trial pe- 
riod between 2.3 to 3.3 d produced similarly disjoint plots, 
suggesting that SU Aur's surface changes on the time scale 
of one stellar rotation. It is likely that we see a mixture of 
line-profile deformations that arise from non-surface events 
and that are due to dynamic short-lived events. 
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4.3 Doppler images 

Despite all these problems, we tried to construct several 
Doppler images over different time spans in the hope of re- 
covering identical features which would allow us to exclude 
shorter lived disturbances. We chose an ephemeris of HJD 
2450409.6. This was picked arbitrarily to correspond to the 
mid-dataset maximum He l absorption. So as to be able 
to differentiate between more persistent structures and very 
short-lived features as well as possible artefacts, we ran test 
reconstructions of data subsets where we excluded certain 
days or also data from individual observatories. 

The reconstruction parameter of the code that we em- 
ploy for the Doppler imaging is the spot filling factor: each 
pixel on the stellar surface is characterized by the amount 
that it is covered in spots (see e.g. Collier Cameron & Un- 
ruh 1994; Barnes ot al. 1998). We assumed the photosplioric 
temperature to be at 5500 K and used the LSD profile of the 
Sun as a lookup table for the immaculate photosphere^. This 
profile is in principle different for each observatory, though 
it turned out that we could use the same template profile at 
the INT and CFHT, provided the lines were only weighted 
by their depths. Whether a depth-only weighting is chosen 
or whether the flux level is also taken into account at the 
CFHT and INT does not have a noticeable effect on the 
resulting images and fits. 

It is not entirelly clear what temperature and line pro- 
file to adopt for the inhomogeneities. On the one hand, the 
relatively constant equivalent width suggests that any ex- 
cess flux, e.g. due to hot spots, is not strongly modulated. 
The erratic changes in the line profile, on the other hand, 
are untypical for cool spots, at least in the framework of 
other active-star observations. We tried three different as- 
sumptions for the inhomogeneities. One where the profile 
changes are due to cool spots, one where they are due to 
hot spots with a "stellar" atmospheric structure, and one 
where the hot spots only emit a continuum black-body spec- 
trum, i.e. where no absorption line is observed for the hot 
structures. 

Depending on which phases and observatories were in- 
cluded in our reconstructions, periods ranging from 2.5 to 
3.1 days yielded the best fits with Doppler images that did 
not show obvious artefacts or very fragmented spot cover- 
age. For the images presented here we chose a period of 
2.7 days as this agrees best with the periods deduced from 
the analysis of the H/3, Na D and He i lines. We further as- 
sumed the inclination angle to be 60° . The inclination angle 
can usually be estimated by determining the best possible 
fits as a function of the inclination angle. In our case, where 
the stellar period is not known exactly, we determined the 
best-fit inclination angle for a range of periods. The angles 
were only weakly dependent on the period and we settled 
on 60° as this yielded good fits for periods around 2.7 d. 
It is also in good agreement with the analysis of IR and 
mm-observations by Akeson et al. (2002) who found an in- 

^ We found that using a different template star (e.g. tfie slightly 
cooler HD 217014) does not change the fits as the equivalent 
width of the profile is adjusted to match the equivalent width 
of SU Aur at each observatory. These adjustments are usually 
very small, large offsets would indicate that the template is not 
suitable. 



clination angle of 621^0* . For our values of the period, rota- 
tional velocity and inclination we calculate a stellar radius 
of approximately 3.6 solar radii (see also Sec. 6). 

4-3.1 Cool spots 

For the cool-spot solutions shown in Fig. 17 we assume a 
spot temperature of 4500 K. The top row shows two images 
produced from data at the CFHT and INT. For the image 
on the right-hand side, the profiles were rescaled to the mean 
equivalent width. Both images are very similar, though the 
rescaled images require a smaller surface coverage to achieve 
an equally good fit. Figs 17(c), (e) and (g) show images 
obtained from observations at the CFHT, INT and OHP. 
While image (c) is for all available data spanning more than 
3 rotation periods, image (e) and (g) are for roughly one 
rotation period centred on HJD 2450407.8 and 2450413.2 
respectively. The contributing observatories and the selected 
time spans for all the images are listed in Tab. 10. Also listed 
arc the values for for the actual fit, the total data set, 
and for a reduced data set where the dimming event and 
the profiles from the BAO which show the lowest S/N have 
been excluded. The fits that result in Fig. 17(d) are shown in 
Fig. 18. The solid lines are the fits and the dashed lines are 
the predicted spectra for the phases that were not considered 
in the fits (due to the dimming event). 

In order to improve the fits to the line profiles, we 
have included differential rotation into our Doppler imag- 
ing code, according to Q{9)/flcq — 1 — Dr sin^(^). In agree- 
ment with Johns-KruU (1996) we find that differential rota- 
tion coefficients of about = —0.1 yield marginally bet- 
ter fits with smaller spot coverage (compare rows 4 and 6 
in Tab. 10). Changes in the differential rotation coefficient 
require adjustments in the value for vsini. For a differen- 
tial rotation coefficient of —0.1 we obtain the best fits for 
wsini = 58 km s~^, while vsini = 59 km s^^ yields the best 
fits if no differential rotation is assumed. For solid-body ro- 
tation, the spot coverage is higher at intermediate to high 
latitudes than for differential rotation. The location of the 
surface features, however, remains the same irrespective of 
whether differential rotation is included or not. This is il- 
lustrated in the Figs 17(d) and (f). Image 17(f) is for a 
differential rotation coefficient of -0.1, image 17(d) for no 
differential rotation. 



4.3.2 Hot spots 

The magnetospheric disk accertion model predicts hot spots 
or bands at high stellar latitudes. If such hot spots are un- 
evenly distributed and are much hotter than the photo- 
sphere, or cover a large fraction of the surface, their ex- 
cess emission should produce a variable dilution of the pho- 
tosplioric line profiles (veiling) in addition to line-shape 
changes. We do not see any consistent trend above the 2.5 %- 
scatter in the equivalent widths of the LSD profiles and pre- 
vious optical veiling determinations have yielded null results 
for SU Aur (see Basri & Batalha (1990); Johns & Basri 
(1995b), though note also that Folha & Emerson (1999) 
found SU Aur to be hightly veiled in the infrared K and 
J band). - 

In a first attempt, we assumed that the spot spectrum 
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Figure 17. Doppler images. The reconstruction parameters and the data sets used for the different images are listed in Table 10, along 
with the spot filling factors and the goodness of fit. The lines underneath the images indicate the phases at which the data used for the 
reconstructions were taken. 
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Figure 18. The LSD profiles along with the fits (solid lines) and the predicted theoretical ijrofilcs (dashed lines) for the different 
observatories. The data are shown as dots along with their error bars. Only the data points overlaid with solid lines were used to obtain 
the Doppler image (shown in Fig. 17d). The dashed lines are the predicted profiles from Fig. 17(d) for those phases that were excluded 
in the fit as they might be affected by a dimming event not related to a change in the surface structure. At the INT, we also excluded 
the very first exposure as it was taken almost a whole rotation period before the next set of observations. The numbers labelling the 
profiles arc the phase of the profile (assuming a 2.7-d period and an ephemeris of HJD 2450409.6) and a rotation count number. The 
observatories were the data, were ta.ken are indicated below each siibfienire. 
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Table 10. List of parameters used for the reconstructions of the Doppler images shown in Fig. 17. The first column lists the identifier. 
Columns 2 to 4 list the goodness-of-fit for different data subsets. Xf gives the goodness-of-fit for the actual fit (using the data as listed), 

while Xt li^t^ ^^'^ goodncss-of-fit between the complete data set and the predicted profiles; X— lists the goodncss-of-fit between the 
predicted profiles and a data set that excludes profiles from BAO and profiles taken between HJD 2450407.75 and 2450409.0 (see also 
dashed lines in Fig. 18). Column 5 shows the spot filling factor. Columns 6 to 10 indicate which time spans of the data have been used in 
the fits for each observatory. The times are listed in modified HJD, i.e. HJD-2450400; a dash indicates that no data from the particular 
observatory was used; "excl" indicates that all profiles were used for a given observatory except for profiles taken between MHJD 7.75 
and 9.0. Image h) is the (thresholded) mean image obtained from images e) and g) and (19) refers to the hot-spot image shown in Fig. 19. 
The last two rows list Xt X?. for unspotted input images, including and excluding differential rotation. 
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Figure 19. Doppler image for continuum hot (6000 K) spots. 
Note that as on previous plots, the photosphere is left white while 
grey indicates the spot filling factor. The total surface covered by 
the hot spots here is 3.3 %. 

would be represented by an ordinary hot-star spectrum. It 

turned out, however, that no satisfactory fits could be ob- 
tained if the convolved hot-spot spectrum showed a signifi- 
cant absorption line. This led us to try fits where the spot 
produces a pure black-body spectrum. Temperatures of 6000 
and 6500 K yielded slightly better fits than the cool-spot as- 
sumption while needing only very small spot-covering frac- 
tions of around 3 %. Fig. 19 shows an example imago ob- 
tained using data from the CFHT and INT. Note that the 
location of the spot features essentially agrees with that 
found from the cool-spot fits (Figs 17(a) and (b)) to the 
same data set. This is not surprising as it is the decreas- 
ing total flux for the cool spot rather than the shape of the 
cool-spot spectrum that acts as the main factor in produc- 
ing the line-profile deformation. Our method does not allow 
us to determine an accurate hot-spot temperature, as the 
size of the deformation is a function of the temperature and 
the surface coverage. We find, however, that the goodness- 
of-fit starts to increase and the reconstructed images start 
to show strongly fragmented spots for temperatures above 



7000 K. For very hot temperatures, the veiling predicted by 
the spots that are invoked to fit the line-profile deformations 
is no longer compatible with the relative constancy of the 
line equivalent width. 



4-3.3 Comparison between different Doppler images 

Whilst the agreement between the different Doppler im- 
ages is not very convincing at first sight, there are features 
that can be found on almost all of the images. To highlight 
those features wc created a "thresholded mean image" of the 
Doppler images displayed in Figs 17(e) and (g). To create 
this image, we only attributed spots to those pixels where 
the spot-filling factor in both contributing images was above 
a certain threshold, in this case greater than 5 %. The fea- 
tures visible on the "thresholded mean image" (sec Fig. 17h) 
are hence features that were present on both input images. 
As image (e) is based on data taken between MHJD 6.4 
to 9.7 whereas image (g) is based on data taken between 
MHJD 11.4 and 15.1, the mean image should show the pix- 
els where structure was present (and possibly survived) after 
more than one stellar rotation had pEissed. 

The mean image in Fig. 17h shows mainly high-latitude 
features. In fact, most images, including Fig. 19 for the hot 
spot, recover two high-latitude spotted areas (between lon- 
gitudes 30° and 112° and from about 270° to 350°) that 
are made up of a couple of stronger spots. Also shown on 
most images are some low- and intermediate-latitude spots 
around approximately 120° and 210° longitude. Note that 
the thresholding biases against low-latitude structure. This 
is because the longitude resolution is latitude dependent and 
much higher at low latitudes. Errors in the period and differ- 
ences in the phase coverage produce small longitude shifts. 
These arc particularly noticeable for narrow elongated low- 
latitude features. The shifts then result in structure being 
suppressed on the thresholded image. 

While we cannot exclude the presence of hot spots with 



stcUar-liko atmospheres, the fits for cool spots or for black- 
body hot spots yield much better goodness-of-fit parameters 
and less fragmented surface images. While we cannot decide 
on the basis of our data whether the spots are hot or cold, or 
whether there is indeed a mixture, the lack of repeatability 
in the photospheric line profiles does favour the hot-spot 
scenario. The lowest high-latitude spot covering fractions 
are seen around longitudes 180° to 210°, i.e. roughly out of 
phase with the maximum He D3 equivalent width and would 
fit in very well with the "eggbcatcr" scenario. 

It is heartening that the mean image is very similar to 
the image reconstructed from all the data save those pro- 
files that might be contaminated by the "dimming event" 
(Fig. 17d). This event is mostly visible for the data at the 
MDO, for which the discrepancy between the predicted line 
profiles and the actual line profiles as shown in Fig. 18 is 
largest. The goodness-of-fit values for the mean image are 
listed in row (h) of Tab. 10. Note that SU Aur's light curve 
shows very erratic behaviour. This may mean that the stel- 
lar surface as we observed it during the MUSICOS campaign 
was not "typical" , if there is indeed a "typical" spot cover- 
age for SU Aur. We therefore do not expect our images to 
show much similarity with previous or future images. 
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LSD profiles of SU Aur, 1996 Nov. 19 
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Figure 20. The Stokes I (bottom) and V (top) least-squares 
deconvolved profiles of SU Aur on Nov. 19. Both profiles have 
been normalised. The Stokes V profile has been shifted by 0.05 
and amplified by a factor of 25 so as to show the structure in the 
line more clearly. 



5 MAGNETIC FIELD MEASUREMENTS 

The exposures (9 groups of 4 subexposures) obtained at the 
CFHT with the MUSICOS spectropolarimeter can be used 
to extract information about the magnetic field structure of 
SU Aur. The first step in this process consists of extracting 
the LSD Stokes V signatures from all available photospheric 
lines. To maximise our chances of detecting the field, we con- 
volved all spectral lines whose relative depth (prior to any 
kind of macroscopic broadening) exceed 10% of the contin- 
uum level in an ATLAS9 G2 subgiant synthetic spectrum 
(Kurucz 1993). As a result, about 3,500 lines were used in 
the analysis. We averaged all spectra obtained on the same 
nights whose rotational phases are within a few % of each 
other thereby decreasing the noise level in the resulting cir- 
cular polarisation profiles further. From the original 9 pro- 
files we thus obtain 5 LSD Stokes V spectra that correspond 
to rotational phases of 0.09, 0.39, 0.83, 0.49 and 0.87, with 
relative noise levels of 0.031%, 0.040%, 0.100%, 0.042% and 
0.037% respectively. The strongest signature we see is ob- 
served at phase 0.09 on Nov. 19 (see Fig. 20). With a peak- 
to-peak amplitude of 0.15% (5a), this signature is associated 
with a detection probability of only 94%. We therefore con- 
sider that this detection is only very marginal, confirming 
that SU Aur is an interesting target for future magnetic field 
investigations with more sensitive spectropolarimeters. 

One can note that the shape of the marginal Zeeman 
signature discussed above, if real, is rather complex with 
at least three polarisation sign switches throughout the line 
profile. This suggests that the associated surface field topol- 
ogy is not simple. It is reminiscent of the polarisation signa- 
tures detected on the weak-line T Tauri stars V410 Tau and 
HDE 283572 with a similar instrument (Donati et al. 1997) , 
that also argue in favour of a rather complex field structure 
such as that reconstructed on the older ZAMS star AB Dor 
(Donati & Cameron 1997). We nevertheless investigated the 
suggestion of Johns-Krull k, Basri (1995b) that the mag- 



netic field of SU Aur is largely dipolar with the magnetic 
axis inclined with respect to the rotation axis. In this model 
accretion occurs mainly at the magnetic poles. Although this 
assumption may sound rather incompatible with the above 
conclusions that the field is complex, it can nonetheless be 
useful to estimate the upper limit to the dipole field strength 
that our data can allow. 

We therefore computed a number of synthetic Stokes V 
spectra for various dipole models, that we compare with our 
9 LSD Stokes V profiles with the help of a simple statis- 
tics'*. While the phase at which the magnetic pole comes 
closest to the line of sight is fixed to that at which the He 1 
line at 587.6 nm exhibits its strongest absorption (i.e. at 
phase 0.0 in the ephemeris mentioned above), the strength 
of the dipole field and the inclination of its axis with respect 
to the rotation axis are considered as free parameters, in 
the range of -1 to 1 kG and 0° to 90° respectively. For all 
these models, we compute the reduced test between the 
synthetic and observed Stokes V profiles, and display it as 
a map in Fig. 21. We can thus conclude that the strength 
of such a dipole field must range between -500 and 100 G 
to remain compatible with our observations (at a 68% con- 
fidence level). Moreover, given the fact that the range of 
reasonable orientations of the magnetic axis with respect to 
the rotation axis is rather 30° to 60° (to ensure that one 
pole is much more visible than the other and that accretion 
is only observed during one half rotation cycle), we finally 
obtain that the maximum dipole field strength ranges be- 
tween about -300 G and 100 G. 

It is difficult to see how such a weak field will be 
able to magnetically confine material accreted at a rate of 
6 X 10~* Mq yr— 1 at a distance of about 2 stellar radii 
(where the local field strength is weaker than 20 G). This 
suggests that while the eggbeater model successfully ex- 

* The results described hold also for the co-added set of 5 LSD 
profiles 
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Reduced map, i=60° 
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Figure 21. Contours of reduced between synthetic and ob- 
served Stolces V profiles. The y-axis sliows the inclination of the 
magnetic to the rotation axis, while the x-axis shows the strength 
of the dipole field. The two darkest gray levels correspond to 
models within the 68% confidence limit, while white corresponds 
to models outside the 99% confidence limit. Field strengths that 
yield acceptable fits for an inclined magnetic axis are rather small 
ranging between about -300 and 100 Gauss. 

plains many of the basic phenomena observed on SU Aur 
and how material is accreted onto its surface, the details of 
the magnetic interation between the central star and circum- 
stellar disc are likely to be more complex. We checked for 
circular polarisation signatures in the He I line at 587.6 nm 
such as those reported by Johns-Krull et al. (1999) on the 
classical T Tauri star BP Tau. No features above the noise 
level are seen in the vicinity of the He I line. 



6 DISCUSSION AND CONCLUSIONS 
6.1 Stellar parameters 

Our analysis confirmed a projected rotational velocity of 
approximately 60 km s~^ for SU Aur, while the rotation 
periods determined from different lines are all below the pre- 
viously determined rotation period of 3 d. Together with our 
estimate of the inclination angle obtained from the Doppler 
imaging procedure this allows us to rc-cvaluate the stellar 
parameters. Unfortunately, the errors on our period deter- 
minations and also on the inclination angle are fairly hefty 
so that we can only give a somewhat large range of probable 
radii. 

For a rotation period of 2.7 d, we find that the 
Doppler imaging process favours inclinations between 50° 
and 70°. The projected rotational velocity of 59±1 km s~^ 
hence translates into an equatorial rotation velocity Vcq = 
68tg° km s"^ Using R = Pv^j2iv, we obtain a stellar ra- 
dius of 3.6^0 3^0. The smallest radius possible for a period 
of 2.7 d is 3.0 Rq. The uncertainty on the stellar period is 
quite large, and periods of 2.5 d and 2.9 d which are also 
plausible would yield smallest radii of 2.8 and 3.3 Rq re- 
spectively. 

Akeson et al. (2002) found that an inner radius of 0.05 to 
0.08 AU for the K-band emission yielded the best model fits. 



This would correspond to about 2.7 to 4.3 stellar radii. This 
is in very good agreement with the usual assumption of the 
magnetospheric accretion model that the inner disk starts at 
the corotation radius. Depending on the exact parameters 
adopted for the stellar mass, radius and rotation period, the 
corotation radius of SU Aur is at about 2.5 to 3 stellar radii. 
Any of the inner disk edges cited above are large enough to 
prevent occultation of the star by the disk. Indeed, for a 
system seen at an inclination of 60° , the disk would have to 
start at 0.7 stellar radii from the surface in order to occult 
the star. 

The derived radii tie in well with the Hipparcos par- 
allax and SU Aur's colour. The Hipparcos catalogue (ESA 
1997) lists a parallax of (6.58±1.92) x 10^'^ arcsccs, a colour 
oi B — V = 0.833 mag and a V-band brightness of 9.23 ± 
0.056 mag. SU Aur's absolute (non-extinction-corrected) V- 
band magnitude is hence 3.3lg;g mag. The bolometric cor- 
rection is independent of the stellar luminosity class and is 
—0.2 mag for the colour term given above (Lang 1980). The 
interstellar extinction has been given as Ay = 0.9 mag by 
Chen et al. (1995) and as Ay = 0.93 ± 0.14 mag by Cohen 
& Kuhi (1979), while the prescription given in Hillenbrand 
(1997) yields Av = 0.46 ± 0.23 mag. Using an extinction 
of 0.9 mag, we obtain a corrected bolometric absolute mag- 
nitude of 2.2^0 ij. For a black body, Mboi (in magnitudes) 
is related to its radius (in Rq) and temperature (see Lang 
1980) through 

Mboi = 42.31 - 5 log - 10 log T. 

Using a temperature of 5700 ± 100 K, we obtain a radius of 
3.2lJ'5_R0. Note that these radii (along with the colour of 
the star) suggest that SU Aur is a subgiant rather than a 
giant. 

The bolometric magnitude of SU Aur translates into 
a luminosity of 10.2^| g Lq. According to the evolutionary 
tracks by D'Antona & Mazzitelli (1994) (using the convec- 
tion description by Canuto & Mazzitelli, 1992), SU Aur's 
mass has to lie between about 2 and 2.5 M© while its age is 
around 3 x 10^ years. This is in reasonable agreement with 
the findings of Dewarf, Guinan & Shaughnessy (1998), who 
determined a mass of 1.9±0.1 M© and an age of 4x 10^ years. 
The large age difference is mainly due to the difiierent ex- 
tinction corrections that were applied. 



6.2 Variability 

SU Aur is certainly a very complex system. While it displays 

some clear periodic features, notably in the Balmer lines, we 
also found short-lived and non-periodical features that make 
the interpretation of the system rather difficult. We find that 
the model by Johns & Basri (1995b) is on the whole success- 
ful in explaining the variability of the Balmer lines, Na D 
and He D3, as well as the correlations between those lines. 
A more in-depth analysis of time lags and cross-correlations 
between the different lines was presented in Oliveira et al. 
(2000b) where some problems of the model along with pos- 
sible extensions were also described. 

During the MUSICOS 1996 campaign, we observed 
transient features in the blue wings of the Ha, H/3 and Na D 
lines. While they appear roughly contemporaneous in all 
lines, their velocity shifts and strengths differ from line to 
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line. Whether the two events arc due to the same mecha- 
nism is not clear as they show different strengths and time 
behaviour (see also Oliveira et al. 2000b). 

Prom the analysis of the Balmer lines, He i (543.6 nm) 
and Na D, wc find strong indications for a stellar rotation 
of the order of 2.8 days, albeit with relatively large uncer- 
tainties, mainly due to the non-sinusoidal nature of the vari- 
ations combined with the relatively short observation time 
span. The picture is less clear for the photospheric lines. The 
LSD profiles do not show any clear periodicities and leav- 
ing the period variable in the Doppler image reconstructions 
yields different optimum periods (spanning 2.4 to 3.2 d) for 
difi^erent data subsets. 

During the campaign, the brightness of SU Aur dropped 
by 0.5 magv. This dimming seems to be coupled with a 
change in the line-profile shape of the photospheric lines. 
The convolved profiles show a redshifted emission features 
for a duration of about 1.2 d, almost half a stellar rotation 
period. The most common explanation for the rapid dirn- 
mings of ETTS is that circumstellar matter obscures the 
stellar disk. This is not normally expected to influence the 
shapes of the photospheric line profiles. We suspect that 
not the complete stellar disk is obscured. There docs not 
seem to be a change in the equivalent widths of the photo- 
spheric lines due to the dimming. The fact that we could also 
not find a wavelength-dependent behaviour of the equivalent 
width suggests that the obscuring matter is grey. 

The dimming is also associated with an increase in the 
equivalent width of the (blue) Ha emission wings, an in- 
crease in the blue emission peak and of the red wing for 
velocities below about 150 km s~^ (see also Fig. 5). If the 
equivalent width of the emission peak is plotted against mag- 
nitude, we find that there is a clear trend of decreasing equiv- 
alent width with increasing brightness. While the brightness 
drops by 0.5 magi^, the equivalent width increases by over a 
factor of two. No such trend is observed for the blue emission 
peak of H/3 even though there is otherwise a good correla- 
tion between the blue wings of Ha and H/3. The fact that no 
line other than Ha shows a correlation with the broadband 
magnitudes suggests that we are not dealing with a contrast 
effect. We would like to note, however, that the S/N is poor 
in the H/3 and photospheric lines, so that available matches 
between well-determined equivalent widths and photometric 
measurements are few. 

While the least-squares-deconvolved profiles contain a 
certain amount of non-repeated and probably also spurious 
structure (see e.g. the line depression at about 25 km s~^ for 
several B AO profiles) , the images reconstructed from various 
data sets all show some consistent high-latitude features as 
shown in Fig. 17h. This is a strong indication that SU Aur is 
indeed covered by a few cool surface spots that last at least 
a few stellar rotations. 

We also want to caution that the surface spots pic- 
tured in Fig. 17h can not be the only features producing 
line deformations. This is illustrated by the only moderate 
improvement of the "goodness-of-fit" compared to what one 
gets from a non-spotted star (see Tab. 10). Note that even 
the stacked residual profiles (leaving apart any period prob- 
lems) show very few events that look like surface spots (sec 
Fig. 15). There are a couple of such events, e.g. between 
about HJD 2450407 and 2450408 and then shortly before 
2450410. But nothing obvious later than that. SU Aur very 



likely presents a mix of periodic (spot or accretion curtain 
induced) variability typical for late-type TTS and the more 
erratic light-level changes typical for UXORs that are prob- 
ably due to variations in the circumstellar dust. Our data 
suggest that some of the UXOR (or typo ill Hcrbst et al. 
1994) variability also affects the shape of the photospheric 
line profiles. 
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